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Abstract 

We study the generation of sharp features in the primordial spectra within the frame¬ 
work of effective field theory of inflation, wherein curvature perturbations are the con¬ 
sequence of the dynamics of a single scalar degree of freedom. We identify two sources 
in the generation of features: rapid variations of the sound speed Cs (at which curvature 
fluctuations propagate) and rapid variations of the expansion rate H during inflation. 
With this in mind, we propose a non-trivial relation linking these two quantities that 
allows us to study the generation of sharp features in realistic scenarios where features 
are the result of the simultaneous occurrence of these two sources. This relation de¬ 
pends on a single parameter with a value determined by the particular model (and its 
numerical input) responsible for the rapidly varying background. As a consequence, we 
find a one-parameter consistency relation between the shape and size of features in the 
bispectrum and features in the power spectrum. To substantiate this result, we discuss 
several examples of models for which this one-parameter relation (between Cs and H) 
holds, including models in which features in the spectra are both sudden and resonant. 



1 Introduction 


Current Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) observations 
favor cosmological initial conditions consistent with a Gaussian distribution of curvature 
perturbations, parametrized by an almost scale invariant power spectrum. In the absence 
of better explanations, these observations are usually considered to support the idea that 
our universe underwent an early stage of accelerated expansion [^, in the form of canonical 
single-field slow-roll inflation [^|^, wherein primordial perturbations evolved adiabatically 
in a quasi-de Sitter background, driven by the evolution of a scalar field gently descending 
the slope of a smooth scalar potential. Such models generically predict an almost Gaussian 
distribution of primordial fluctuations with a slightly red tilted power spectrum in 

addition to a small value for the tensor to scalar ratio, in agreement with the latest CMB 
constraints [^[^. 


It is possible, however, that primordial fluctuations had their origin under less placid 
conditions, with their dynamics characterized by time-scales shorter than the inverse of the 
quasi-de Sitter expansion rate H, breaking the simple statistical properties offered by single- 
field slow-roll scenarios. If present, such time-scales are expected to be associated to the 
existence of massive degrees of freedom (with masses larger than H) either affecting the 
inflationary background or interacting with curvature perturbations as they exit the horizon 
during inflation [^. The landmark of models incorporating time-scales shorter than 
is the emergence of features in the primordial spectra, with their size and shape charac¬ 
terizing the physics associated to the high energy degrees of freedom that participated in 
their production [9 . Consequently, the existence of features in both the primordial power 


spectrum and bispectrum could represent a unique chance of acquiring a better understanding 
of inflation by giving us insights into mass scales —additional to H — characterizing the 


fundamental theory responsible for them 17 


There are well known examples of inflationary models allowing for features in the pri¬ 
mordial spectra. One obvious class corresponds to canonical single-field models in which the 
inflaton field traverses a sudden change in the slope of its potential, for example in the form 
of a step or a bump, resulting in a brief interruption of the slow-roll dynamics. This type 
of dynamics leads to a sudden variation of the expansion rate, without necessarily implying 
large deviations from the quasi-de Sitter background as the leading geometrical configura¬ 
tion during inflation. These sudden variations of the background are inevitably felt by the 
primordial fluctuations while they exit the horizon, leaving a set of features imprinted in 
the power spectrum and bispectrum (18-30 . The same mechanism leads to the production 


of features in P{X, (f) models if the (^-dependent part of the Lagrangian has a feature that 
interrupts the slow-roll evolution of the field as it traverses the target space of the theory. 
In this case, not only the expansion rate will display a sudden variation in time, but also 
the sound speed c<j at which curvature perturbations propagate (typically encountered in 
non-canonical realizations of inflation such as P(X, i;/))-models), leading to a richer variety 
of features imprinted in the primordial spectra as compared to those produced in canonical 
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models of inflation. 

Another class of models allowing for features in the primordial spectra are multi-field 
models of inflation]^ in which the inflaton background trajectory meanders the landscape 


offered by the multi-field target space 37-41 , enhancing the interaction between curvature 


perturbations and fields orthogonal to the trajectory each time the trajectory undergoes a 
bend [8l|liy42ll45| . The resulting interactions have the chance of generating localized features 


in the spectra if the bend-rate is larger than the Hubble expansion rate H during inflation 46 
In this case there are two main regimes leading to features: 


I. Features generated by multiple degrees of freedom: If the time variation of the bend-rate 
of the trajectory is comparable to —or larger than— the mass of the fields orthogo¬ 
nal to the trajectory, these become excited, generating features in the spectra with 


a wavelength characterized by the mass of the excited degrees of freedom 11,47 48 


These degrees of freedoms may be heavy or light, depending on whether their masses 
are larger or smaller than the inflationary expansion rate H. 

II. Features generated by a single degree of freedom: On the other hand, if the time varia¬ 
tion of the bend-rate is smaller than the mass of the helds orthogonal to the trajectory, 
then the only dynamically relevant degree of freedom consists of curvature perturba¬ 
tions [^[^. It is then possible to study these features using an effective single-field 
theory describing the evolution of curvature perturbations, in which the background 
is parametrized by the sound speed Cg of curvature perturbations and the expansion 
rate H (4^59| . In this case, the features in the primordial spectra are found to be the 
consequence of the time variation of Cg and H resulting from the bend in multi-field 
space iBSIll- 

While in the second regime the background is of the multi-field type, the dynamics associated 
to the fluctuations maps into an effective single-field theory description. This is possible 
because the heavy degrees of freedom may decouple from the light curvature perturbations 
even if the bend rate is much larger than the expansion rate 11 . Nonetheless, the multi¬ 


field dynamics associated to the heavy degrees of freedom continue to be present through 
the background, leading to the emergence of nontrivial couplings in the effective field theory 
for the fluctuations, such as the sound speed Cg describing the propagation of curvature 
perturbations. 

Despite of the complexity involved in the emergence of features, the study of models allow¬ 
ing for them may be systematically organized —and substantially simplified— by adopting 


the effective held theory (EFT) framework (M,63 in which one disregards the details of 


the background dynamics to focus, instead, on the dynamics of curvature perturbations and 
any other degrees of freedom relevant around horizon crossing. To be more specihc, the 


*Yet another mechanism leading to the emergence of features —also involving additional degrees of 
freedom— is offered by particle production during inflation 36 . We will not be concerned about this 
type of models in this article. 
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EFT approach allows one to study the generation of features in the primordial spectra by 
parametrizing any nontrivial dynamics containing tiny time-scales into background couplings, 
such as the sound speed c* and the expansion rate H appearing in the EFT Lagrangian for 
the fluctuations. This is precisely the case of both P{X, (p) models and multi-field models 
with heavy fields (second regime already mentioned), where the emergence of features in the 
primordial spectra may be traced back to sudden time variations of the sound speed Cg and 
the expansion rate H (and subsequent time derivatives) appearing in the effective Lagrangian 
describing the dynamics of curvature perturbations. 


1.1 Correlating features in the primordial spectra 


The purpose of this article is to study the emergence of features in the primordial spectra by 
employing the effective field theory of inflation formalism. Adopting this framework will allow 
us to uncover model independent relations between features appearing in different n-point 
correlation functions of primordial curvature perturbations, in particular between features 
in the power spectrum and the bispectrum. This will turn out to be possible because the 
EFT action of the curvature perturbation depends only on a reduced number of background 
parameters present at every order of its non-linear perturbative expansion. In this way, by 
relating the time-variation of these parameters, we will be able to deduce a set of consistency 
relations correlating features in the inflationary primordial spectra. 

For instance, it is already well known how to express the bispectrum as a function of the 
power spectrum in the particular case of features due to small sudden changes of the sound 
speed Cg of curvature perturbations away from the value = 1 (see also ref. [^). The 
existence of such a relation is allowed by the possibility of expressing the sound speed as a 
function of the power spectrum in the following way 



dk AV 

~kl^ 




( 1 . 1 ) 


where r is the usual conformal time, and the difference AV = V — Vo parametrizes features 
in the power spectrum V{k) against its featureless counterpart V(){k) (which has a slight red 
tilt). This equation results from analyzing the relation determining the power spectrum V{k) 
as a function of Cs(r) and inverting it. On the other hand, given that the combination 1 — 
appears as a coupling in the cubic part of the EFT action for the fluctuations (generating 
non-Gaussianity), having Cg as a function of the power spectrum permits one to compute the 
bispectrum in terms of the power spectrum. Specifically, one finds that the /nl parameteij^ 
has the general dependence on the power spectrum given by 


fNLikl,k2,k3) = 


;^^dlnk^ Vo '^^Ulnk Vo ^ Vo 


k —{fci+fc2H“^3)/2 


( 1 . 2 ) 


where fi = fi{ki,k 2 ,k 3 ) are known smooth scale-independent functions of the triangle con- 
hguration ki -|- k 2 -|- ks = 0, determined by the fact that features come from sound speed 


^The /nl -parameter is defined in eq. (2.301. 
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variations 64 . Equation (1.2) is found to reduce to the well known Maldacena’s consistency 


relation [66| - |68| in the squeezed limit, where one of the momenta is much smaller than the 
other two {i.e. <C ki,k 2 )- Moreover, in the particular case of sharp features, where the 
sound speed is subject to rapid variations with a characteristic time-scale much smaller than 
the features become dominated by the second order In/c-derivative, giving us back the 
following general expression for sharp features 


/nl ^ P 


(P AV 
dink‘s Vo 


(k) 


. k=(k\+k2+k3)/2 


(1.3) 


where j3 = /3(fci, k 2 , kp) is a smooth function, independent of the overall scale, but dependent 


on the shape of interest and the source of the feature (/2 in the parametrization of eq. (1.2)). 


It is possible to generalize eq. (O) to the case in which features are generated by both, 
time variations in the sound speed and time variations in the expansion rate H. This gene¬ 


ralization was derived in ref. 69 , and is found to be given by 


Ifl _ PY'" _L ^ = A 
8^ 2r2 r4 ni 


C+oo 




dk k-^^{k)e-‘^^^^, 
) '0 


(1.4) 


where 6h — and t] = e/He (with e = —H/H^ the usual slow roll parameter and ' 

standing for derivatives with respect to conformal time r). However, because both 1 — and 
T] are couplings appearing in the cubic action, it is clear that under more general backgrounds 
the bispectrum cannot be expressed uniquely in terms of the features appearing in the power 
spectrum, inevitably introducing a degeneracy in the parameter space determining the source 
of features]^ Despite of this, eq. (1.4) may be used to express the bispectrum as a function of 
the power spectrum in two limiting cases: when features are generated by a variation of the 
sound speed, and when they are generated by a variation of the expansion rate. To be more 
specific, in the case of sharp features, if features are generated by a variation of the sound 


speed one finds that the bispectrum has the form shown in eq. (1.3) with (3 given by 


/3s = I 


kik2ko fef -|- ^2 + ^3 


Qk\ + k^ + kl {ki + k2 + k^Y ’ 


(1.5) 


whereas if features are generated by a variation of the expansion rate, then (3 acquires the 
form given by: 

kik2ko 


Pe = 


( 1 . 6 ) 


12 k\ + kl + kl' 

These two results for (3 may in principle be tested by searches of features in the CMB non- 
Gaussian angular spectrum. 

where a general expression for the bispectrum in terms of 


Mhis point was specially emphasized in ref. 
the power spectrum was obtained. 


70 
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1.2 Article’s main idea &: results 


Equation (1.3), together with the functions (3^ and /?£ of eqs. (1.5) and (1.6), parametrizes 


features in the two limiting cases in which these are generated by a variation of the sound 
speed and the expansion rate respectively. From these results, it is reasonable to expect that 


eq. (1.3) continues to be valid under more general circumstances in which these two sources 


occur simultaneously, as expected in more realistic models. As we will show in this article, 
this is indeed the case: there exists a natural one-parameter realization for the function j3 


in eq. (1.3) such that /3s and j3^ of eqs. (1.5) and (1.6) constitute the desired limiting cases. 


This one-parameter realization of /3 is a consequence of a relation between rj, parameterizing 
the evolution of the expansion rate, and the sound speed c* given by 


ado 

V = Vo + 


(1.7) 


where a is a constant specified by the model underlying the emergence of features, and rjQ 
represents the slowly varying part of rj. Although it is clear that this non-trivial relation is not 
guaranteed to be valid in general, we will argue that it remains a fairly good approximation for 
a large variety of models admitting sudden time variations of their slowly varying background 
quantities, including P{X, 4>) and multi-field models with heavy fields. For the purposes of 


the present work, admitting the validity of eq. (1.7) enables us to interpolate between the 


two limiting cases shown in eqs. (1.5) and (1.6), obtaining a much more general relation for 
/3 given as: 


/3a = 


1 


hhks 


kf + k2 + A:| 
121 + akf + k^ + " {ki + k 2 + k^Y 


a -|- 2 


( 1 . 8 ) 


This relation gives us back eq. (1.5) in the specific case of a = 0 and eq. (1.6) in the 
case |a| —>■ -|-oo. Equation (1.8) constitutes the main result of this work. It consists of a 


single parameter consistency relation between the power spectrum and the bispectrum valid 
for triangle conhgurations away from the squeezed limit, where one recovers Maldacena’s 
consistency relation. It offers a concrete —and fairly model independent— parametrization 
of features to be tested in future primordial spectra reconstructions using cosmic microwave 


background and/or large scale structure observations 79-88 . 


1.3 Outline 

The present article is organized as follows: We will start in Section where we propose 


the use of eq. (1.7) as a valid prescription relating the time dependence of both the sound 


speed and the inflationary expansion rate. There we show that by admitting this relation 


one obtains eq. (1.8), linking together the features in the primordial bispectrum and the 


primordial power spectrum. In Section we justify the use of eq. (1.7) and examine its 
validity within different held theoretical realizations of inhation involving the presence of a 
sound speed. To be more precise, we numerically solve inhationary backgrounds with sudden 
features in various classes of P{X, cj)) and multi-held models of inhation. Then, in Section 
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we show that our proposed relation (1.7) also serves to parametrize resonant features. Finally, 
in Section we provide our concluding remarks on the findings of this work. 


2 A general parametrization of features 


The appearance of sharp features in the primordial spectra may be traced back to the time- 
variation of background quantities at a rate much larger than those characteristic of slow- 
roll inflation (73| - |78| . To analyze the effect of these rapid variations on the dynamics of 
fluctuations, we adopt the EFT of inflation perspective. This framework asserts that the 
action describing the dynamics of curvature perturbations may be parametrized with the 
help of a limited number of background parameters, determining non-trivial relations between 
coefficients appearing at different orders in the theory. In comoving gauge, the action for 
primordial curvature perturbations, up to cubic order, may be written in the following way 

S = + (2.1) 


where the quadratic part of the action, corresponds to (written in units such that 

= 1 ) 


5 ( 2 )= d^xah 


^7t2 _ i_(v7^)2 


whereas the cubic part, 5(^), is given by 


89 


5(3) = 


1 


d^x a^e 

c 

l-c2 


^ [3(c^ - 1) + € - ??] 


(1 -c2) + 7? + e- 


^ 7?3 _L ^2v72tt. ^ ^V72 






eo'^ 


( 2 . 2 ) 

^ TZ{vnf 

(2.3) 


2cs 

H^sJ 

f <55(2) 


<57^ 


where x is given by the constraint equation V2x = o^eTZ/c^. In these expressions, the 
expansion rate H and the slow roll parameters e and rj parametrize the evolution of the scale 
factor a{t) as 

rr ® ^ ^ A\ 

^“a’ "'“(He’ 

where the dot represents a derivative with respect to cosmic time t. A relevant background 
quantity parametrizing coefficients appearing in both quadratic and cubic parts of the action, 
corresponds to the sound speed Cs{t) which determines the speed at which long wavelength 


modes propagate. On the other hand, the parameter A in eq. (2.1) parametrizes the strength 


of the operator TZ^, and it is usually found to depend on Cs according to a relation determined 
by the specific model in question. Finally, the quantity / multiplying the linear classical 
equation of motion 53^'^'^/6TZ is a given quadratic function of TZ, whose specihc form will turn 
out to be irrelevant for the present discussion. 

This action may be simplified substantially by making a few reasonable assumptions 
about the type of features we want to study. First of all, we will consider variations of the 
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sound speed Cg such that it stays close to the value Cg = 1. That is, we will assume that 


0 = 1 - < 1 , 


(2.5) 


at all times. We hnd this assumption reasonable, because a suppressed sound speed Cg is 


known to produce a sizable level of equilateral non-Gaussianity 90091 , so far undetected 71 


An immediate consequence of this choice relates to the size of the coefficient in front of the 
operator appearing in eq. (2.3). Indeed, from an effective field theory point of view, it is 


possible to show that this coefficient is in fact naturally of order 9^ (see for instance, ref. 92 


1-c" 


2A 


(l-c; 


, 2^2 


9 ^ 


( 2 . 6 ) 


Another important assumption that we will adopt is that the time variation of the background 
does not break the quasi de Sitter phase characterizing inflation during the horizon crossing 
of the modes responsible for the primordial spectra accessible to observations. In other words, 
we will assume that the evolution of H is such that 


e <C 1, 


(2.7) 


at all times. In practice, because e is restricted to be positive, eq. (2.7) implies that e will 
have the form 

e = cq + Ae, (2-8) 

where cq corresponds to the slowly varying part of e, characterizing the average slow-roll 
behavior of inflation, and Ae is the part containing rapid departures from cq. In addition to 


condition (2.7), we will assume that Ae is such that its variation is small relative to eo- In 
other words: 


Ae <C eo. 


(2.9) 


This condition means that the background of our interests contains rapid variations, but will 
continue to stay close to the quasi-de Sitter stage parametrized by the featureless value eo- 
The previous splitting of e implies that r] may be written in terms of Ae in the following way 


1] = Vo + Ar/, Arj = 


Ae 


( 2 . 10 ) 


where rjo = eo/Heo represents the slowly varying part of r]. Finally, we will assume that the 
time variation of both Ae and Cg are characterized by a time scale much smaller than 
(which is what we mean by rapid variations). Consequently, the features we are interested 
in are characterized by the following hierarchies: 


1 

H 


d9 

dt 


»», i 


dAe 


dt 


> Ae. 


( 2 . 11 ) 


It may be appreciated that condition (2.11) implies that Ai], defined in (2.10), may be of 


order 1 (or even larger) without implying a breaking of condition (2.7). More generally, we 
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assume that any background quantity A parametrizing rapid variations (such as 0, Ae and 
Ary) satisfies: 


1 

H 


dA 

dt 


»l^l- 


( 2 . 12 ) 


It is also useful to cast this hierarchy in terms of both conformal time r and e-folds N. Let 
us recall that conformal time is dehned to satisfy the relation dr = dt/a, whereas N satisfies 
dN = Hdt. Then, the previous hierarchy may be written in the following two alternative 
ways, that will turn out to be useful in the next discussions: 


dA. 

dr 


» 1^1, 


dA 


dN 


»l^l- 


(2.13) 


To continue, all of the previous assumptions may be put together to organize and simplify 
the form of our initial action given in eq. (2.1). In particular, the cubic part of the action is 
found to be dominated by the following terms: 


c(3) _ 

•^int ~ 



{39 + r])nn‘^ + 


4(r0'-ry)7^(V7^)2|. 

J 


(2.14) 


This result allows us to study the effects of features in the bispectrum using the in-in for¬ 
malism of cosmological perturbation theory. Details of how this is done may be found in 
ref. |69| . Here instead we quote the main results of Appendix in order to derive various 
relevant expressions to obtain the primordial spectra. 


2.1 Features in the power spectrum 


Let us first consider the effects of a rapid variation of the background on the power spectrum. 
The dimensionless power spectrum V(k) is usually defined in terms of the 2-point correlation 
function of curvature perturbations in Fourier space (evaluated at the end of inflation) as: 

(7^k7^k') = (27r)35(k + k')^r(k). (2.15) 

The computation of (T^kl^k') may be performed using the in-in formalism of perturbation 
theory discussed in Appendix [A| In particular, one finds that the power spectrum P(k) has 
the following form (see Appendix [A]) : 

Po = ^(^) . (2.10) 


Here Vq corresponds to the zeroth order featureless power spectrum, which coincides with the 
conventional prediction offered by canonical single field slow-roll inflation (n^ is the spectral 
index parametrizing the small scale dependence of the power spectrum, and fc* is a pivot 
scale fixed by mode that exited the horizon when the background was characterized by an 
expansion rate Hq). On the other hand, AV is the part containing the features, and is given 


by 69 


AV 

-^(k) = k 

Vo ^ ^ 


!\r 

' —CO 


-9 + 


6h , 2(5_h- 




1 ddn 


k‘^T‘^ k'^A k^A dr 


sin(2A:T), 


(2.17) 
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where r represents conformal time. In addition, 6h is a function containing information 
about the time varying background, which in the case of rapid variations is found to be given 


by 69 


A 1 ' 

= -^rri . 


(2.18) 


We may quickly verify the validity of eq. (2.17) for the particular case in which 9 = 6q 1 
is a constant and 6h = 0. In this case the integral of eq. (2.17) may be solved by dropping 
the highly oscillatory part of the lower limit r —?• — oo, giving us back 


AV 


00 

2 


i-ci 


1 - Cs. 


(2.19) 


This result, together with eq. (2.16) gives us back the result V = Vojcs = HQ/Sn'^eoCs, which 
is the well known expression valid for the power spectrum for curvature perturbations with 
a constant sound speed. 

To continue, eq. (2.17) may be simplified by assuming that both 0 —)• 0 and d/f —)• 0 as 
T —)■ 0“, which simply means that small variations of the background happening at the end 
of inflation are meaningless for the purpose of computing features in the long wavelength 
power spectrum, describing modes that exited the horizon long before inflation finishes. This 
assumption allows one to perform several integrations by part and write: 


,AV 


1 


k^ — {k) = -- 
7^0 


dr 


- 0 "" 


W _ o . 


2r2 


sin(2/cr) 


( 2 . 20 ) 


Let us emphasize here that up to this point, we have made no assumptions about the sharp¬ 
ness of features. This new expression may be Fourier-inverted by conveniently extending the 
time domain of the functions 6 and 6h to the whole range r G (—oo, -|-oo), and demanding 
them to be odd under r —)• —r. The result is found to be given by eq. (1.4) already shown 
in the introduction. However, if we further use the fact that 6h satisfies the hierarchy of 
eq. (2.13), we are allowed to drop the lower order time-derivatives from the left hand side of 
eq. (1.4), from where we hnally obtain: 
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T]_ 

4r 


d 




dk k^‘^{k)e-^^^^. 
Vo ^ ^ 


( 2 . 21 ) 


This relation gives us the time dependence of a specihc combination of background quantities 
in terms of the features appearing in the power spectrum. Before moving on to consider 
features in the bispectrum, it is important to notice that the hierarchy of eq. (2.13) together 
with eq. (2.21) implies a hierarchy between fc-derivatives for the power spectrum, which takes 
the form: 


d^ AV 


<1 


AV 

dink‘s Vo 


dln/c Vo 


Vo 


( 2 . 22 ) 

In other words, rapid variations of the background (characterized by a timescale much smaller 
than H~^) imply sharp features in the power spectrum. We will use this important result in 
the next discussions, and verify its validity in Section where we examine several concrete 
examples of models with features. 
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2.2 Features in the bispectrum 

We may now examine the appearance of features in the bispectrum. Let us recall that 
the bispectrum B may be defined through the 3-point correlation function of curvature 
perturbations, in Fourier space, at the end of inflation as 

(T^kiTtka'^ka) = (2vr)^(5(ki-|- k2 -I- k3)i?(ki, k2, ks), (2.23) 


where, because of the delta function (i(ki -|- k 2 -|- ka), the three momenta are restricted to 
add up to zero. By computing the three point correlation function using the in-in formalism 
one finally obtains 

B = Bo + AB, (2.24) 

where Bq represents the standard featureless, slow-roll suppressed bispectrum [^. On the 
other hand, AB is the part containing features, and is found to be given by (see Appendix 

0 


AB(ki,k2,k3) = 


27r^V^ 


dr i X 


J- 

(30 + rj) [-(A: 2 A: 3 )^ - (kskif - {kik2f + ikik2k3T{k2k3 -|- k^ki + kik2)] 

7] — tO' 


+ 


2r2 


[kl + k2 + fcf] (1 - ikiT)(l - ik2T)(l - ik^r) 


(2.25) 


with the shorthand K = ki + k 2 + k^. This relation is equivalent to that of eq. (2.17) and 
allows us to compute the shape and size of features in the bispectrum out from the varying 
background quantities 6 and rj. 


2.3 A proposal to relate background quantities 


Crucially, both AV and AB depend on 6 and rj simultaneously. This implies that AV and 
AB are in principle not tied together, and may contain scale dependent features which do 
not correlate with one another. Nevertheless, we may break the degeneracy in the para¬ 
meter space of these observables by assuming a dynamical relation between the background 
quantities 9 and rj. In this respect, our proposal consists of assuming the following relation 

ry 

ri = m- (2.26) 


where ryo = ^o/H eq, and a is a slowly varying dimensionless function that may be considered 
to be a constant for all practical purposes. We will justify the use of this relation in the next 
section, where we will verify its validity within a wide range of models, including P(X,(p) 
and multi-field models of inflation. In the meantime, notice that this relation, together with 
eq. (2.20), implies that the power spectrum is now given by: 


AV 


(k) 


— (l + a)k f dr 9 {t) sm(2kT). 

J —OO 


(2.27) 


10 








On the other hand, using eq. (2.21) allows us to determine the form of both 6 = 1 — Cg and 


At/ = r/ — ryo in terms of the features appearing in the power spectrum as: 


9 = 


1 


AV 


Arj = 


1 + a vri 
1 


a 


OO 

H-OO 


k 


{k)e 


1 + a TTZ 


dk 


Vq 
d AV 
dkl^ 


—2ikT 

(2.28) 

(A:)") 

(2.29) 


It may be seen that in the limit a —)• 0 one recovers the case in which the features are 
exclusively due to variations of the sound speed, and in the limit |q:| —)> +oo one reobtains 
the case in which features are due to variations of the expansion rate. One may worry about 
the limit a —)• —1, which makes eqs. (2.28) and (2.29) diverge. In that limit, subleading 


terms that were neglected at the left hand side of eq. (2.21) take over, and modify the final 


forms of eqs. (2.28) and (2.29). In the present analysis we omit this specihc situation and 


leave the related problem of finding such expressions open (we further comment on this issue 


in Section 3.3). 


2.4 A general correlation between featnres 

Now that we have both 9 and r] expressed in terms of the piece AV determined by features. 


we may use eq. (A.34) to reach a general expression determining AB as a function of AV. 


To obtain it, it is enough to replace eqs. (2.28) and (2.29) back into eq. (A.34). To write the 


result, it is convenient to introduce the standard /NL-parameter as a function of ki, k2, and 
ks, defined as: 

10 kik2k3 {kik2k^)‘^ 


/nl = 


3 kl + kl + kl {2 tiYVI 


-AB. 


(2.30) 


Then, by keeping the leading contributions to AB, which according to eq. (2.22) consists of a 


term proportional to second order derivatives of AV in terms of In k (and a term proportional 
to a first order derivative that dominates in the squeezed limit) we finally find 


/nl 


/3o 




dink‘s Vo 




d AV 


(k) 


12 dink Vo J 


(2.31) 


where /3a = /3a{ki, /c 2 , ko) is a scale-independent function given by: 


13a = 


kik2ko 


a -|- 2- 


kf k2 3- 


12 1 3- a kf 3- k^ 3- k^ (ki 3- k2 3- k^Y 


(2.32) 


This final form of /? interpolates between the results and (|1.6[) discussed in the intro 


duction, found in the two separate cases in ref. 69 . In particular, it is worth highlighting 


that in the equilateral configuration /ci ~ /c 2 — A: 3 , /3 has the form 

5 1“ 


= 


36 1 + a 


2 

a+- 


(2.33) 
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while in the folded configuration ki = k 2 = k^/l, f3 becomes: 

1 1 


/5^foid) ^ 


3 

“+4 


(2.34) 


12 1 + a 

Finally, it is worth noticing that in the squeezed limit (fei = ^ 2 , k^ —)• 0) we recover the well 
known consistency relation (see also ref. [^), as the coefficient of the second derivative of 
the power spectrum disappears 

d AV 


f(sq) _ 

Jnl - 12 


dink Vo 


(k) 


(2.35) 


independently of the value of a. As already mentioned in the introduction, eq. (2.32) const! 


tutes our most important result. It provides a unique relation determining the shape and size 
of features in the bispectrum in terms of those appearing in the power spectrum with only 
one parameter to adjust. As we shall see in the next section, the parameter a is determined 
by the model responsible for features, and in principle may be constrained by observations 
aiming to characterize the scale dependence of the the primordial spectra. At any rate, it 


is important to recall that eq. (2.32) is strictly valid away from the value a = — 1, and that 


it assumes that features were generated by rapidly varying background quantities satisfying 


eq. (2.13). 


2.5 On the general strnctnre of the correlated spectra 


The results of this section ratify the general dependence anticipated by eq. (1.2) of the 


introduction. A striking aspect of this relation is the fact that /nl does not depend on third 


order derivatives of AV/Vo or even higher. In fact, (1.2) constitutes a truncated version of 
a more general relation between /nl and derivatives of AV/Vo that may be expressed as: 


n=oo 


/NL(fei, A:2, ko) = y /, 


n=0 


dJ^ AV 
dlnk^ Vo 


(k). 


(2.36) 


However, it is possible to verify that fn with n > 2 are slow-roll suppressed. To appreciate 
this, first notice that In/c-derivatives acting on AV/Vq are the consequence of terms contai¬ 


ning powers of conformal time r appearing inside the integral of (A.34). These in turn, come 


from the time dependence of the curvature perturbation wave function TZk{T) given in (A.32), 
valid in the de Sitter limit e —)> 0 (see the discussion of the Appendix . Nevertheless, this 
specific form of TZk{T) receives corrections of order e that introduce additional powers of r 


inside (A.34), finally leading to (2.36) but with coefficients fn (with n > 2) suppressed by 


slow-roll parameters. This observation gives us an upper limit on the sharpness of features 
that our method is allowed to study. That is, in order for the truncation to remain valid, 
features in the power spectrum must be such that: 

-/2 A-n 

(2.37) 

In other words, our approach is useful to analyze sharp features, but their sharpness is limited 
by the size of the slow-roll parameter e. 


d^ AV 


<1 

iM 

d^ In k Vo 


(V In k Vo 
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3 Sound speed and expansion rate during inflation 


Our main result (2.32), deduced in the previous section, depends crucially on the validity 


of eq. (2.26) relating the sound speed Cg and the slow-roll parameter r]. In this section we 


provide arguments to support this conjecture and discuss a few examples for which it is 
found to be valid. Let us start by emphasizing that, a priori, there are no reasons to expect a 
simple relation determining the value of r] in terms of the sound speed Cg since, in principle, 
the background dynamics has enough degrees of freedom to allow for situations in which rj 
and Cg evolve independently. Nevertheless, in the class of models that we are attempting 
to describe, features appear as the consequence of small deviations from the quasi-de Sitter 
background driven by inflation, characterized by constant values of e and c^: 


e = eo, 


2 2 
= Co- 


(3.1) 


In these models, e parametrizes the quasi de Sitter geometry whereas Cg parametrizes the 
non-trivial kinematical properties of the fluid that is causing the quasi-de Sitter dynamics in 
the first place. It is therefore reasonable to envisage that a small change in e away from cq 
will come together with a small compensating change to the sound speed away from the 
value Cq. In other words, whatever may be the cause of the rapid background variations, it 
should induce simultaneous variations of e and c^. As long as these variations are small, we 
expect that 

AeocAc^, Ae = e-eo<Ceo, Acf = cf - Cg <C Cg, (3.2) 


where eg and cg are the slowly varying values of e and Cg, that may be taken as constant values 
for all practical purposes. The idea behind this relation is simple: each time a constant value 
of e is reached (that is, each time ^Ae = 0) a new quasi-de Sitter geometry is achieved, 
which should be characterized by a new constant value of the sound speed Cg. This precisely 
requires that both quantities vary in synchrony. The proportionality relation Ae oc Ac^ 
allows us to deduce a more useful relation between r] and Cg. Indeed, by using eq. (3.2) back 


in eq. (2.10) where we defined Arj, we obtain the non-trivial relation 


V = Vo 


a dAc^ 
2 dN ’ 


(3.3) 


where a represents some slowly varying background quantity determined by the specific 
model causing the rapid variation. Then, in the particular case in which the sound speed 
varies away from the background value cg = 1, we finally obtain 


V = Vo ^ 


a dO 
2dN' 


(3.4) 


which is our desired result. In what follows we check the validity of these results by analyzing 
several concrete examples. 
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3.1 Features in P(X, 0)-models 


For our first two examples we consider a class of /c-essence models with a Lagrangian 
of the form C = P(X, (/>), where P{X,(j)) is a given function of X = —{d4>)‘^/2 and cj). The 
homogeneous background equations describing the evolution of inflation are given by the 
Friedman equation and the equation of motion for the scalar field cj), which, in terms of 
number of e-folds N, are found to be 


3P2 -E = 0, 


Cp(j) 

dN^ 


I /q 2 ^ ^ 4 ^ I 


P2p 


P(A — 0, 


X 


(3.5) 

(3.6) 


where now X = <^^/2 (due to the homogeneity of the background), E = 2XPx — P is the 
energy density of the inflaton field, and Px = dxP, etc. In these models, the 

slow roll parameter e and the sound speed Cg are respectively given by: 


XPx 2 ^ Px 
p 2 ’ Px + 2XPxx 


(3.7) 


Both of these quantities stay almost constant in a quasi-de Sitter regime. Thus, the origin of 
features may only appear as the consequence of a non-trivial dependence of these quantities 
on the rolling inflaton field (j), producing displacements of e and away from their quasi-de 
Sitter fiducial values eo and Cq . These displacements will happen not only due to the explicit 
appearance of (p in eq. (3.7), but also on the time variations of X responding to this non¬ 
trivial dependence. However, as we shall see, as long as these displacements are small, they 
will happen in synchrony, in such a way that quasi-de Sitter is recovered in terms of both 
parameters at the same time. In what follows we examine a few examples for which this 
happens and corroborate the validity of eqs. (3.2)-(3.4). 


3.1.1 Example 1: Localized features in the potential 


Let us start by considering the specific case in which a rapid variation of the background is 
caused by a localized feature in the potential. To be concrete, we choose a P{X, (^)-theory 
of the form 


P{X, c/)) = X + AX^-V{(l)), 


(3.8) 


where A parametrizes a non-trivial contribution to the kinetic term, quadratic in X, and 
V{(p) corresponds to a chaotic potential with a small feature on top of it in the following 
manner 




(3.9) 


where m is a mass scale parametrizing the featureless chaotic potential, and f{4>) is a function 
parametrizing the localized feature. To discuss this class of models let us consider two 
particular choices for f{4>). For our first choice, model-(a), we consider a small step feature 
of the form 


f{a){(p) = Htanh 


V2Ap\ ’ 


(3.10) 
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where (^o gives us the position of the step, A4> its width, and B the size of the jump produced 
by the step. For our second choice, model-(6), we consider a small Gaussian bump in the 
potential of the form 


= Bexp 


(3.11) 


2 A (()2 

where, again, cpQ represents the position of the step, Acp its width, and B the amplitude 
of the bump. These potentials have been studied in the work 1^. To be consistent with 


the hierarchy of eq. (2.13) we assume that both features are “sharp” in the sense that their 


effects on the solutions 4>{N) and H(N) take place within an e-fold of inflation. In practice, 
this means that A(j) should be such that Afp <C \d(p/dN\. In this way, even if we choose the 
amplitude B of the features such that they have a tiny effect on e and c^, the sharpness due 
to the small value of A(p may induce a dramatic effect on rj and d^O. More to the point, let 
us consider the following values (in Planck units) for the parameters in both models 


A = 10^°, <po = 14.19, B = 5x Acp = 0.002, m = 5.975 x 10 


-6 


(3.12) 


These parameter-values imply a background sound speed Cq ~ 0.8. It is important to assert 
that the formulas of the previous sections —involving the primordial spectra— are strictly 
valid for the case Cq = 1, and therefore we are not allowed to use this example to infer a 
relation between the bispectrum and the power spectrum as discussed in Section In spite 


of this fact, this example is still useful to reinforce our conhdence on the validity of eqs. (3.2) 


and (3.3). We have chosen the value of 4>o in such a way that it produces a feature about 60 


e-folds before the end of inflation, roughly implying the generation of features in the power 
spectrum within the window of scales relevant for CMB observations. 

Figure summarizes the main numerical results obtained by solving the equations of 
motion (|3.5l) and (3.6) for the two models in question. The left hand side panels show the 


results pertinent to model-(a), whereas the right hand side panels show the results for model- 
(6). In particular, the top panels show the functions Ae = e — eo and Ac^ = — Cq as a 

function of N. On the other hand, the bottom panels compare Ar] and as a function 

of N for the values 


a 


(a) 


15, 


a 


(b) 


15.6, 


(3.13) 


for models (a) and (b) respectively. We see that the plots ratify, in a rather eloquent way, the 


validity of eqs. (3.2) and (3.3). The features on the potential with the parameters values of 
eq. (3.12) imply the rapid variation of background quantities within a window of AN ~ 0.075 
e-folds. In both examples the variations experienced by Ae = e — eo and Ac^ = — Cg 

remain small, but the variations of Arj and djyd are found to be large. Moreover, the fact 
that a^a) 3.nd are larger than 1 implies that the features are mostly the result of the 
rapid variation of the expansion rate. We have tried values for the parameters different from 


®In all our numerical examples, we define the end of inflation by e = 0.1, and choose onr parameters such 
that "Po = 2.43 • 10“® 60 e-folds before the end of inflation. The featnre takes place inside the CMB-window: 
in all plots the point A'’ = 0 corresponds to 55 e-folds before the end of inflation. 
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Figure 1: The figure shows the main numerical results for model-(a) (left-panels) and model-(6) (right- 
panels). The top panels show the functions Ae = e — eo and Ac^ = — Cq as a function of N. The bottom 

panels show At/ and ac?jv0/2 as a function of N. Notice that to satisfy eq. (3.31, model-(a) requires a ~ 15 
whereas model-(6) requires a ~ 15.6. 


those of eq. ( |3.12 ) and have found similar results, except for the cases in which conditions 
Ae <C eo and | Ac^| 


<C Cq are violated. 


3.1.2 Example 2: DBI inflation with featnres 

For our second example we will stay within the realm of P{X, (/))-models, and consider a DBI 
theory 9^ 97 parametrized by the following Lagrangianj^ 


P{X,Tp) = f-\^) [l - Vl - 2 /(</))a1 -E(</.). 


(3.14) 


Here, /(</>) corresponds to the warping factor describing a system of D3 branes on a warped 
background accomplished in certain string compactifications. /(</>) couples together the ra¬ 
pidity (j) and the vacuum expectation value cf), allowing for departures from = 1 depending 
on the shape of /(</>). Observe that in order to ensure that we stay in the physical regime 

^An alternative example, which renders similar results to those offered by this example, consists of 
P{X, 4>) = X + g{4>)X^ — V{<}>), where g{4>) accomplishes the same role of /(<(>). 
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1 — Cg > 0 at all times, we require /((/>) > 0 everywhere. For definiteness, let us consider two 
choices for f{4>) characterized for being suppressed everywhere except for a region of finite 
support. Our first choice, model-(c), consists of a Gaussian bump of the form 


/(c)(</>) = ^exp 


i4> - 

2A(t>‘^ \ ’ 


(3.15) 


where and A(/> are the position and width of the bump, and B is the amplitude modulating 
the coupling between (p and X. For our second choice, model-((i), let us consider the following 
form for consisting of a Gaussian bump multiplied by {(p ~ 


f(d0) = B 


A(?i2 


exp 


{(p - (po)^ ' 
2 A (/>2 


(3.16) 


Here again, cpQ and Acp are the position and width of the bump, which this time consists 
of two peaks around po- These two alternatives for /(</>) induce a brief exchange between 
both kinetic and potential energy when the vacuum expectation value of p hits the value po. 
Moreover, the non-trivial dependence of the Lagrangian on X implies a sound speed different 
from unity during the duration of this exchange. To examine the effects of these sudden 
variations of the background, let us consider again a chaotic potential of the form 


ViP) = 



(3.17) 


able to realize inflation within canonical single field slow roll inflation, with 60 e-folds of 
inflation, and adjust the values of po, Ap, and B to generate variations such that Ae <C eo 
and 1 — Cg <C 1. We have chosen (for both models) 


Po = 15.387, B = 5x 10^ m = 5.875 x 10"®. (3.18) 


Moreover, in model-(c) we set Ap = 0.005 while in model-(d) we have Ap = 0.004. Figure]^ 
shows the main numerical results for our two models. The top panels show the background 
solutions for both Ae = e — eo and 1 — Cg as a function of e-folds N. It may be seen that they 
happen in synchrony, as already anticipated by our general discussion linking both quantities. 
The bottom panels show both Arj = rj — r]o and a /2 as functions of N. We find that 

there is a good agreement between both quantities for the values 


«(c) - -0.53, a(rf) ~ -0.52. 


(3.19) 


The resulting plots ratify the validity of our ansatz (3.4). We have checked that this relation 
is ensured as long as the conditions Ae <C eo and 1 — Cg <C 1 are satisfied. As soon as these 
conditions are broken, one starts to observe a disagreement between both quantities. 


3.2 Features in multi-field infiation 

Let us now take a look into the appearance of features within the framework of multi-field 
inflation. Concretely, we consider a multi-field model with a set of two scalar fields (/>“, 
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Figure 2: The figures show some relevant results concerning models (c) and (d) (left and right panels 
respectively). The top panels show 1 — cl and Ae = e — eo as a function of N. The bottom panels show 
Arj = Tj — rja and adNd/2 as a function of N. 


a = 1,2, with a non-trivial sigma model metric 'jab, and a Lagrangian given by 

£ = -l^abdrd<p^-V{<p), 


(3.20) 


where is the scalar potential responsible for producing inflation. Given that we are 

focusing our interest on the case of multi-field models with only two helds, it is convenient to 
introduce a pair of unit vectors T“ and tangent and normal to the inflationary trajectory 
(f)°- = dehned as 


Ta^^dr 

00 dt 

Na = Vdet 


(3.21) 

(3.22) 


where Cab is the Levi-Civita symbol in two dimensions (notice that these definitions keep the 
orientation of the two vectors hxed). In addition, we have defined 


00 = 


7fec0'’0''- 


(3.23) 
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As usual, the metric '^ab and its inverse 7 “^ are used to lower and rise indices on tensor living 
in the tangent space of the multi-field manifold parametrizing the model. The vectors T“ 
and N°‘ allow us to introduce the local rate of turn of the inflationary trajectory n defined 
to satisfjfJ 

DT°- 

-^ = -nN^, (3.24) 

where /dt = T°‘ + is the covariant time derivative of T“ with being the 

Christoffel symbols computed with the help of 7 afe. The quantity O is nothing but the 
angular velocity of the tangent vector T“ as the inflationary trajectory bends. To continue, 
the equations of motion dictating the evolution of the fields deduced from the action 


(3.20), are given by: 


D(j)° 

dt 


+ 3 /^ 0 “ + V^ = 0, 
3H^ = \4 + V{cP). 


(3.25) 

(3.26) 


The two equations of motion for the scalar fields (3.25) may be projected along T“ and 
respectively, giving us back 


(^0 + 3H(^q + — 0, 

4>o 


(3.27) 

(3.28) 


where = r“14 and Vn = N°'Va. On the other hand, by differentiating eq. (3.26) with 


respect to time, eq. (3.27) allows us to show that: 


_ iio_ 

2H^' 


(3.29) 


A detailed analysis of the dynamics of the two scalar fluctuations in this class of theory shows 
that perturbations may be decomposed into curvature perturbations and a massive degree of 


freedom, associated to the fluctuations perpendicular to the inflationary trajectory [37-39 
The mass m of the massive degree of freedom is found to be 


= Vnn + 


+ 3f2^ 


(3.30) 


where Vnn ^ N<^N^(Vab - ^ab^c)) and M is the Ricci scalar computed out of ^ah- Then, it 
turns out that if m? ^ it is possible to deduce a low energy effective theory for curvature 


perturbations alone, of the form ( 2 . 1 ), where the sound speed is given by [ 8 ]: 


c7' = 1 + 


4122 


m? — 412^ 


(3.31) 


I Notice that this quantity has been denoted 9 in refs, [s 11 studying the effect of turns on the dynamics 
of curvature perturbations. Here we opt to use Q, in order to avoid any confusion with 6 = 1 — c7 
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In other words, the sound speed will always decrease as soon as the inflationary trajectory 
suffers a turn. 

Now that we have a concrete relation giving us the sound speed in terms of quantities 
parametrizing the inflationary trajectory, we can understand the validity of relation (3.4) 
under a new scope. First, it is important to appreciate that departures of the sound speed 
from the value Cs = 1 are exclusively due to turns of the inflationary trajectory n 7 ^ 0. In 
addition, eq (3.28) tells us that each time there is a turn, the trajectory is pushed against the 
wall of the potential away from the minimum Vjv = 0 (which, from a Newtonian point of view, 
it is a consequence of the angular momentum accompanying the turn). Then, each time there 
is a turn, the inflationary trajectory will tend to climb the wall of the potential, implying 
an increase in the value of the potential, and a consequential reduction of the kinetic energy 
X = <^o/2- This has the cost of producing a small reduction of the value of e, as implied by 
eq. (3.29). Therefore, a turn simultaneously produces both a reduction of the speed of sound 
Cg < 1 and a decrease of the value of e, finally leading to eq. (3.3). 


3.2.1 Example 3: Localized turns in a two-field target space 


We now proceed to examine a few concrete numerical examples of inflationary trajectories 
with turns, corroborating eq. (3.3). Let us start by choosing the following notation for the 
two scalar fields: 

r = {x,^). (3.32) 


We will consider a model where turns occur as a consequence of a feature appearing in the 
sigma model metric given Xab- To simplify matters, let us parametrize the metric in the 
following wa}** 


1 f{x) 
fix) i + /^(x) 


(3.33) 


where / is a given function of the first field y to be specified shortly. Notice that this metric 
has a unit determinant by construction. In addition, we choose a potential with the following 
generic form 

V(x, V’) = v(x) + (3.34) 

where v(x) is a single field potential constructed to realize slow-roll inflation in the y- 
direction. On the other hand, the term with represents a mass term forcing ?/) to stabilize 
at the value if = 0. Notice that if fix) = 0 the two fields would remain decoupled, and in¬ 
flation would happen along the y-direction, as the y-field descends the slope of the potential 
vix), with if stabilized at the value if = 0. More generally, if / is a constant, then one can 
always redefine the fields and deduce that the inflationary trajectory consists of a straight 
path in the multi-field target space. To make things more concrete, let us consider a simple 

**This type of model was already considered in ref. |ll| to study features in the power spectrum as the 
consequence of turns in multi-field models. 


20 









chaotic potential of the form 


(3.35) 


v{x) = 


m 


-X 


and choose initial conditions in such a way that inflation lasts for about 60 e-folds. Next, let 
us consider two different choices for the profile of the function f{x): characterized for being 
constant everywhere except for a region of finite support. For our first choice, model-(e), let 
us consider the following function: 


/(e) (x) = ^ (^1 -Ftanh 


ix - Xo) 

I V2Ax J 


(3.36) 


This represents a step function, with a width Ay centered at xo that changes from / = 0 
to f = B. It is easy to see that this transition implies a single turn with non-vanishing rate 
of turn n inducing a non-trivial departure of away from unity. For our second choice, 
model-(/), we consider a Gaussian function of the form 


/(/)(x) = ^exp 


(X - Xo)" 

2Ax2 


(3.37) 


This form of / induces two consecutive turns in opposite directions. As a consequence, we 
obtain (for the same set of parameters) a non-trivial departure of away from unity in the 
form of two consecutive bumps. To solve the inflationary dynamics of these two models we 
consider the following values 


B = 2.5, 


= 10 


-2 


m = 5.89 • 10“®, Ax = 0.005, xo = 15.37. (3.38) 


Figure shows the main numerical results for these two models. The top panels show the 
background solutions for both Ae = e — eo and 1 — as a function of e-folds N. It may be 
seen that they happen in synchrony, as already anticipated by our general discussion linking 
both quantities. On the other hand, the bottom panels show both Ar] = r] — rjo and a 8^0/2 
as functions of N. We find that there is a very good agreement between both quantities for 
the following values of a: 

a(e) ^ —0.53, a(^f) ~ —0.52. (3.39) 


Yet again, these numerical results emphasize the accuracy of our ansatz (3.4) to describe 
inflationary backgrounds with rapid deviations from quasi de Sitter. 


3.3 Accuracy of the power spectrum 

Our previous examples (2 and 3) consisted of models that required values of a that were 
worryingly close to the value a = — 1. Indeed, let us recall that our expression for the power 


spectrum given by eq. (2.27) stops being valid for values of a close to —1, where other sub¬ 


leading terms start to dominate the final form of the power spectrum. Nevertheless, we find 


that for the aforementioned examples the power spectrum computed with (2.29) constitutes 


a very good approximation of the power spectrum computed directly by numerically solving 
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Figure 3; The figures show various relevant results related to the models e and / (left and right panels 
respectively). The top panels show 1 — cl and Ae = e — eo as a function of N. The bottom panels show 
Arj = Tj — rja and adNO/2 as a function of N. 


the mode functions exactly. Figure compares the power spectrum obtained by solving the 
mode functions of curvature perturbations exactly (continuous red curves) and our analytical 
approximation (dashed blue curves) using eq. (2.27), for models (c), (d), (e) and (/) respec¬ 
tively. These plots confirm that the approximation offered by eq. (2.27) gives an accurate 
representation of features in the power spectrum, in spite of the value a ~ —0.5. 

The main reason for the good agreement between both power spectra is the sharpness of 
the features. Indeed, let us recall that eq. (2.27) is derived from (2.20) by keeping the leading 
terms according to the hierarchy (2.13). This tells us that the power spectrum of eq. (2.27) 
is accurate as long as 


1 cx 
AN 


> 1 , 


(3.40) 


where AN is the e-fold scale characterizing the rapidly varying background. If this condition 
is not fulfilled (for instance, because a ~ —1), then eq. (2.20) tells us that the leading order 
term becomes: 

AV 5a 6”' 


Of course, this form of the power spectrum modifies our main result (1.8), a situation that 
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Figure 4: The figure shows a comparison between the power spectrum obtained by solving the mode functions 
of curvature perturbations exactly (continuous red curves) and our analytical approximation (dashed blue 


curves) using eq. (2.271. 


will not be examined in the present article. 


4 Resonant features 


In this section we use our formalism to consider the case of resonant features. By resonant 
features we mean features in the power spectrum and bispectrum induced by an inflatio¬ 
nary background with a certain periodic time-dependence 98-101 . We wish to show that 


our ansatz (3.4) continues to constitute a good characterization of the time-dependence of 


resonant departures from quasi-de Sitter, and therefore our relation (2.32) offers a concrete 
parametrization for the search of resonant features in the primordial spectra. Before consi¬ 


dering a concrete example, let us first study the implication of eq. (3.4) analytically in the 


case of resonant features. To do so, let us consider a background characterized for having a 
varying sound speed with a harmonic dependence on the e-fold variable N: 


eiN) = eo [l + sm{uoiN-No))] 


(4.1) 
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Here wq is the frequency of the oscillations (in e-fold units) and Oq corresponds to the average 
value of 0{N). Equation (3.4) then implies that r]{N) is given by 


a 


r]{N) = r]o + -u!o9ocos{ujo{N - No)). 


(4.2) 


We may use these relations to obtain an expression describing the features in the power 
spectrum. Using N = Nq — ln(r/ro), which is valid for a constant expansion rate H, we find 


that (2.20) becomes 


(fe) — -I- (1 -I- a)9ok 


dr sin(u;o ln(r/ro)) sin(2A:r), 


(4.3) 


where the first term, 9o/2, comes from the non-oscillatory piece of eq. (4.1). In fact, our 
discussion around eq. ( |2.19 ) tells us that we may absorb this piece in the definition of the 
featureless power spectrum Vo —t Vs = Volco, where cq = \/l — 9o. Then, integrating 


eq. (4.3) we finally obtain 

AV 


— (A:) = (1 + a)ypocosh(7ra;o/2) sin (^a;oln(2fe|ro|) + 


(4.4) 


where po and (po are real functions of loo defined through the algebraic equation poV''^° = 
r(l —iwo)) where T correspond to the usual T-function. This result may be more conveniently 
summarized as 


AV 

{k) = A sin [cuo In 2k + ip], 

S 


(4.5) 


where p = po + wq in I'^ol- This relation tells us that the power spectrum has a logarithmic 


dependence on the scale k. This result may be plugged back in eq. (2.31) to obtain an 
expression for the bispectrum containing resonant features 


/nl ^ -^a{ki, k2, kz)A(xil sin [wq ln(/i;i + /c2 + ko) + p] , 


(4.6) 


where (5a is given by eq. (2.32). This result may be compared with previous parametrizations 


used to search resonant features in the CMB data [102 , such as the one used by the Planck 
satellite to test resonant features in the non-Gaussian bispectrum [^, with the generic form 
given by 

Bres{ki,k 2 ,ko) = g{ki,k 2 , A: 3 ) sin [C In (/ci + k 2 + ko) + ((] (4.7) 

where g{ki, k 2 , k^) is a function of the shape of the triangle conhguration. In our case, the a 
parameter plays an important role in determining the shape of ^(/ci, A; 2 , /cs). 


4.1 Example 4: Resonant features in multi-field inflation 

Let us now consider an example where the background offers a periodic dependence on time. 


To be concrete, we consider again the multi-field model of Section 3.2.1, where the features 
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Figure 5: The left panel of the figure compares the background quantities »7 — ? 7 o and adj^Q 12. The right 
panel compares the exact power spectrum computed numerically, with the analytical expression given by 


eq. (4.41. 


are introduced via the function /(x) appearing in the ct- model metric of eq. (3.33). In this 
occasion we shall consider the following form for the function f{x) 

.{x - XoY 


fix) = B cos { 2tt 


Ax 


(4.8) 


and keep the chaotic potential u(x) = nnrclr/2 for the inflationary sector x within the multi- 


held potential (3.34). The choice (4.8) ensures that the background quantities appearing in 


the EFT of curvature perturbations will inherit a periodic dependence on time. In particular, 
one hnds that both Cs and e evolve periodically —and in synchrony— just as anticipated by 
eq. (3.2), and one verihes the validity of eq. (3.4). To offer an example, we have used the 


following values for the parameters of the present model: 

B = 1.5, = 10■^ Ax = 0.05, m = 5.975 x 10"®. 


(4.9) 


These values ensure that the departures of Cg and e are small from their respective featureless 
background values 1 and eq- Figure compares both rj — rjo and adjqO/2 for the value 
a = —0.53, which is found to give the best match between the two functions. Moreover, we 


hnd that (4.1) constitutes a good fit of 6{N) with Oq and uiq given by: 


6»o ~ 0.016, Wo ^ 35.4. 


(4.10) 


With these values, it turns out that our expression for the resonant power spectrum of 


eq. (4.4) offers an accurate representation of the effects of the oscillatory time dependence 


of the background. The right panel of Figure compares the power spectrum of the model 


computed numerically with the expression given by eq. (4.4), showing a very good agreement 
between both functions. 


It is important to mention that we have also verified that eq. (4.4) offers an accurate 


expression for classes of models other than multi-field. Indeed, we have analyzed numerically 
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the example of Section 3.1.2 for which we have chosen the warp factor f{4>) to have the same 
form of fix) in eq. (4.8), and have found similar results to those presented here. 


5 Discussion and conclusions 


While it is hard to conceive that observations could ever reveal the nature of the physics 
responsible for features, it is clear that, if they are present at all, both the power spectrum 
and bispectrum will contain complementary information about the source that generated 
them. Indeed, in this work we have seen that, under certain reasonable assumptions, it is 
possible to correlate features appearing in the bispectrum with those appearing in the power 
spectrum. Specifically, we have studied the generation of sharp features by adopting the 
perspective of effective field theory of inflation, where the primordial spectra of curvature 
perturbations are the consequence of the dynamics of a single scalar degree of freedom. In 
this framework, sharp features are the consequence of the simultaneous variations of the 
sound speed Cg and the expansion rate, parametrized by e. 


The novelty of this work is that we have shown that two different kind of features, the 
ones caused by a perturbation in e and the ones caused by a perturbation in Cg , are in fact not 
separate at all. After all, they follow from the same perturbed path that the inflaton takes 
through field space. In this way, in order to correlate the primordial power spectrum with the 
bispectrum, we had to uncover a universal relation linking Cg and e, given by eq. (1.7). This 
relation —valid for small but rapid variations of Cg and e— led to our main result, given by 
eq. (2.31) supplemented with eq (2.32). To show the validity of eq. (1.7), we offered several 
numerical examples of models allowing for sharp features. 


The aim of this work therefore is to break the degeneracy between features in the power 
spectrum induced by perturbations in e, and features in the power spectrum induced by 
perturbations in Cg. Looking at the power spectrum alone, they cannot be told apart. With 
this work in hand one can, however, link the features in the power spectrum to those appearing 
in the bispectrum. One can extract a value for a, which indicates exactly to what extent the 
feature is caused by a perturbation in e, and to what extent by a perturbation in Cg. 


Our results leave several questions unanswered that deserve to be addressed. Let us 
mention a few possible directions that we find interesting: 


1. One may consider the general problem of using the present methods to study the 
appearance of features in higher order n-point correlation functions. The present results 
suggest that features in higher n-point correlation functions may be expressed in terms 
of features in the power spectrum. More concretely, it would already be interesting to 
learn how the tri-spectrum depends on features appearing in the power spectrum. 


2. As emphasized in Section 3.3 our main result (2.31) stops being valid for values of 
a violating condition (3.40), with a the proportionality parameter a between r] and 
dO/dN . Thus, a direction for further research could be to identify and analyze models 
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in which a is sufficiently close to —1. In this case, we need to reconsider the relation 
between the perturbations in the power spectrum and rj and 0, as terms which are 
subleading for a ^ 1 become dominant. 


3. In this work, a is just a constant that has been adjusted to satisfy eq. (3.3) in each 


model. Nevertheless, it is clear that a is a model dependent quantity. Thus, it might 
be possible to deduce analytical expressions for a in terms of parameters characterizing 
the models hereby examined {i.e. P{X,(f>) and multi-field models of inflation). 

4. We have not considered models in which the slowly varying part of the sound speed 


Co remains far from 1. Our first numerical example (Section 3.1.1) shows that a back¬ 


ground sound speed different from 1 does not spoil the relation between ry and dO/dN. 


However, in this situation there are new operators in the cubic action (2.14) for TZ 


inducing features in the bispectrum, not taken into account in the final result (2.31). 


5. To correlate the features in the spectra we have assumed a standard single field effective 


theory (2.2) describing the dynamics of curvature perturbations, in which the sound 


speed is the main parameter describing departures from canonical single field inflation. 
However, as already noticed in [62] , there are more general classes of effective field 
theories where the dispersion relation of curvature perturbations contains non-trivial 
dependences on the scale 55,58 . Therefore, it would be interesting to analyze the 


way in which the primordial spectra are correlated within other classes of effective field 
theories. 

6. We have focussed our efforts exclusively on the case of sharp features, where the hierar¬ 


chy of eq. (2.12) was assumed to characterize every rapidly varying background quantity. 


We see no obstructions towards obtaining a more general expression correlating features 
in the bispectrum with those of the power spectrum in the case where the variation of 
background quantities is not so restricted (as in the case of ref. (^, for the particular 
case of features emerging exclusively from a varying sound speed.). 


We leave these interesting challenges as open problems. 

To conclude, if some of the not-so-significant-yet features that Planck is seeing now turn 
out to be real, our work provides a framework to analyze them. If not, it can be equally 
useful in order to confirm the paradigm of plain, featureless single field slow roll inflation. 
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A Perturbation theory for features 


In this appendix we summarize the in-in formalism of perturbation theory in quasi-de Sit¬ 
ter space-times [^, 103| that we use to compute the power spectrum and bispectrum with 
features. More details may be found in ref. |69| and references therein. 


A.l The in-in formalism for featnres 


As a starting point, let us consider the action (2.1), together with eqs. (2.2) and (2.3), 


describing the dynamics of curvature perturbations TZ in comoving gauge. Then, it is useful 
to introduce the canonically normalized curvature perturbation u{t, x) through the following 
field reparametrization: 

u = zTZ, z = V^—. (A.l) 

Cs 


Doing this change of variable allows us to split the quadratic part of action (2.1) as = 


, where the zeroth order part Sq^’ and the interacting parts , are respectively 

given by: 


:( 2 ) 


^(2) 


( 2 ) 


0(2) _ 1 
*^0-2 

0 ( 2 ) _ 1 
int 2 


d^x dr 
d^x dr 


.12 


- {Vuf + 
9{t){Vu)‘^ -I- -^5{t)u 


where 6{t) is a rapidly varying background quantity, given by: 


S = 6h-t 6' + -t 


A 1 ' 

= -^rrj . 


(A.2) 

(A.3) 

(A.4) 


( 2 ) 

Notice that the zeroth order piece, Sq , corresponds to the conventional quadratic action 
describing a scalar perturbation in a de Sitter space-time sliced with the help of cosmological 
coordinates. On the other hand, contains coefficients induced by the rapidly varying 
parts of the background. Of course, we also have the cubic contribution to the total action. 


which, to leading order in the rapid background quantities, is given by (2.14). This means 


that, with the purpose of employing the in-in formalism, we may finally split the entire theory 
as 


5 = + 5int, 


where S'int is the interacting part of the action, up to cubic order, and given by 


o _ c(2) ci(3) 

‘Jint — -I- , 


(A.5) 


(A.6) 


where S'); is given by eq. (|A.3l) and S'bV is given by eq. (|2.14l). Equation (|A.2l) informs us 


:(3) 


that the n-field corresponds to the canonically normalized adiabatic curvature perturbation 
with a canonical momentum given by vr = dCjdu' which, to leading order, reduces to vr = u'. 
Both u and vr satisfy the equal-time canonical commutation relation 


[u(x, r), 7r(y, r)] = i5{x - y). 


(A.7) 
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where (5(x — y) is the Dirac delta function. The field u{x,t) and its canonical momentum 
7r(x,r), satisfying this commutation relation, may be written in terms of the interaction- 
picture field evolved in time with the help of the propagator U{t) as: 

u(x, r) = 17^'(r)u7(x,r)[/(r), (A.8) 

7r(x, r) = C/^(r) tt/(x, r) 17(t). (A.9) 


Here the interaction-picture field ui{x, r) has the form of a free field, which may be expanded 
in Fourier modes as 

1 


u/(x, r) = 


(2vr) 


d-^k 


akUk{T)e"'^'^ + alul{T)e 


zkx 


(A.IO) 


where oj^ and are creation and annihilation operators satisfying the usual commutation 
relations for particle states in a Fock space: 

[ok, flp] = (27r)^(5(k - p). (A.ll) 

On the other hand, Mfc(r) represents the normalized solution to the linear equation of motion 
in momentum space derived from (A.2): 


- n]uk = 0. 


(A.12) 


The solution to this equation is obtained by selecting the Bunch-Davis vacuum, and is given 
by: 


Ukir) = 


1 




1 - 7 ^ ] g-ifcr 

kr , 


Returning to (A.8) and (A.9), the propagator is given by 

1 


U{t) = T exp l—i dT'Hi{T') , 
I J —00 + 


(A.13) 


(A.14) 


where T stands for the standard time ordering symbol and oo+ = (l-l-ze)oo is the prescription 
isolating the in-vacuum in the infinite past. Furthermore, Hj is the interaction picture 
Hamiltonian given by 

Hi{t) = C/t(r)i/intD(T), (A. 15) 


where Ffint is the interaction Hamiltonian derived from (A.6). In particular, it is possible to 
show that the quadratic contribution to the interaction-picture Hamiltonian has the form: 


= I d^x 


e{T){Vuif + 


(A.16) 


(3) 

On the other hand, the cubic contribution to the interaction picture hamiltonian, Hj , is 
found to be given by 


= rup; f d^x a^e 


{w + v)nin‘j + - ri)ni{vnif 


(A.17) 
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( o\ 

Notice that we have opted to express Hj in terms of TZj = —ujHqt/ \/2eo rather than uj. 
These expressions allow us to compute the two- and three-point correlation functions for u 
perturbatively. To first order in the rapid background quantities, one hnds: 

(u(x,r)u(y,T)) = (0|n/(x,r)n/(y,r)|0)-Fz f (iT'(0| [i7/(r'), u/(x, T')tt/(y, r')]|0), (A.18) 

J —OO 

(u(x,r)u(y,T)n(z,r)) = i f dr'(0| [iff T')tt/(y, r')zz/(z, t')]|0). (A.19) 

J — OO 

These expressions conform the basis to compute the power spectrum and bispectrum with 
features, which is done in what follows. 


A.2 Primordial spectra with features 


In order to define the power spectrum, and later on the bispectrum, we write the adiabatic 
perturbations in Fourier space by introducing the mode function 7^k('7‘) as: 


7^(x,r) = 


1 


(2vr) 


d^kiZ\^{T)e 


ikx 


(A.20) 


Then, the dimensionless power spectrum r), evaluated at a given time r, is defined via 

the equation: 

(7tk(r)7tp(r)) = {2TTf6{k + p)‘^Vn{k,T). (A.21) 


Putting together (A.l) and (A.21) we deduce an expression for the power spectrum V'ji{k, r) 

lly no 

(M(x,r)u(0,r))e“**"'^, 


in terms of the two-point correlation function of the canonically normalized tt-held 

k^ 


'Pn{k,T) = 




(A.22) 


where stands for f d^x. Notice that one of the fields has been conveniently evaluated at 
y = 0, allowed by the homogeneity and isotropy of the background. Similarly, the bispectrum 
i? 7 ^(ki, k 2 , ks, r) is conventionally defined as: 


(7^kl(r)7^k2(r)7^k3(T)) = (27r)^(5(kik 2 k3)B7^(ki, k 2 , ks, r). 


(A.23) 


This expression may be inverted to give the bispectrum in terms of the three-point correlation 
function of the n-field as 


57e(ki,k2,k3,T) = \ 


(n(x,r)t6(y,r)zz(0,r))e 


—iki-x—zk2-y 


(A.24) 


X Jy 


where k 3 = —k 2 — ki. Just as we did with eq. (A.22), we have evaluated one of the fields at 
the comoving coordinate z = 0. Finally, let us recall that our goal is to compute correlation 
functions in the long wavelength limit |rA:| <C 1, which gives us the spectra at the end of 
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inflation. Thus, the power spectrum and bispectrum we are interested in correspond to the 
following formal limits: 


'P'R.ik) = lim Vn{k,T), 

hKo 

57e(ki,k2,k3) = lim S 7 e(ki,k 2 ,ks,r). 

hKo 


(A.25) 

(A.26) 


Given that TZ is constant after horizon crossing, these expressions give us the initial condi¬ 
tions for adiabatic perturbations outside the horizon for the hot Big-Bang era, in terms of 
background parameters at horizon crossing time, during inflation. 


To compute the primordial power spectrum with features we just need to plug (A. 18) 
into (A.22). This leads to the following form of the power spectrum 


V = Vo + AV, 

where Vq corresponds to the featureless power spectrum given by: 


87r2eo' 


(A.27) 


(A.28) 


In fact, if we had considered slow-roll corrections to the zeroth order mode, we would have 
obtained the more accurate expression (used in the present work) 

ris — l 

, (A.29) 


k 


^ 8^0 Vfe 

where n* is the spectral index parametrizing the small scale dependence of the power spec¬ 
trum, and A;* is a pivot scale fixed by mode that exited the horizon when the background 
was characterized by an expansion rate Hq. To continue with eq. (A.27), AV is the part 
containing the features, found to be given by 

AV 


— (fc) = lim i{2clk^T'^) dr e [77/(r'), u/(x, r)u/(0, r)] |0) 

To fc|r|-:>0 J_oo Jx 


= k 


dr 


1 ddn 

k'^r'^ ' k'^T^ k^T^ dr 




sin(2/cT) 


(A.30) 


where stands for f d?x. Finally, we may compute the bispectrum in a similar way by 
plugging (A. 19) into (|A.24[). This step gives us 


2eo 


AB(kl,k 2 ,k 3 ) = T^7^l(0)7^2(0)7^3(0) / dr 
^^0 J —c>o 

t 9' — rj 


W + ri 


\JZi{t)'R! 2 {t)TZ' 2 ,{t) + sym\* + c.c. 


[k 2 • ks 7^l(r)7^2(T)7^3(r) -h sym]* -h c.c.^, (A.31) 


where TZi{T) = 7^(kj,r) is the wave function for comoving curvature perturbations obtained 

(A.32) 


from (A. 13): 


Ho 
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The previous relation for the bispectrum may be further simplihed. On the one hand, notice 
that 

' ^ ^ ^ (A.33) 


1 


ki • k2 + k2 • ks + ka • ki = + A:| + A:|), 


which is valid as long as ki + k 2 + ka = 0. Then, by using eq. (A.32) and assuming that the 
functions 9 and r] are odd functions of conformal time r. Then we get 


Afi(ki,k2,ka) = 




dr i X 


{kihks)^ J-c 

[39 + rj) [-(/c 2 fca)^ - {k^^kif - {kik 2 f + ikik 2 k 3 T{k 2 k 3 + k^ki + kik 2 )] 
V-T-9' r;2 


+ 


2r2 


[kf + k2 + fc|] (1 — ikiT){l — ik2T){l — iksT) 


(A.34) 


where K = ki + k 2 + k^, and where we have identihed Vq = //Q/Svr^eo. These are the results 
needed to compute the final expressions linking features in the bispectrum with those of the 
power spectrum. 
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